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We compute the temperature dependence of the antiferromagnetic order parameter and the gap in 
the two dimensional Hubbard model at and close to half filling. Our approach is based on truncations 
of an exact functional renormalization group equation. The explicit use of composite bosonic degrees 
of freedom permits a direct investigation of the ordered low temperature phase. We show that the 
Mermin- Wagner theorem is not practically applicable for the spontaneous breaking of the continuous 
spin symmetry in the antiferromagnetic state. The critical behavior is dominated by the fluctuations 
of composite Goldstone bosons. 
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The low temperature phase of the Hubbard model 
is highly interesting for a variety of phenomena observed 
in strongly correlated electron systems. This includes the 
antiferromagnetic state for zero or small doping, on which 
we concentrate in this note, as well as possible super- 
conducting states for large doping. There has been a 
long debate about the effective theory for the spin wave 
fluctuations 0, ^'^'i correct implementation of 
mean field theory^ the problem being to reconcile 
the S0(3)-spin symmetry with the Hartree-Fock results. 
More recently, functional renormalization group studies 
have made progress in the understanding of instabilities 
in various channels and the corresponding phase diagram 
0) S 0) ^ ) but have been unable to explore the ordered 
phases for low temperature. (For other approaches see 
[llLll3| .l We concentrate in this note on the two dimen- 
sional Hubbard model with next neighbor interactions. 
There a further problem arises for the understanding of 
the antiferromagnetic state - a theorem by Mermin and 
Wagner |13] forbids the spontaneous breaking of a contin- 
uous symmetry at nonvanishing temperature. This seems 
to contradict an antiferromagnetic order since the latter 
breaks spontaneously the continuous spin symmetry. 

In this letter we present a combined study of 
Schwinger-Dyson (or Hartree-Fock) equations [13 . [l5l |. 
extended mean field methods using a general- 

ized Hubbard-Stratonovich transformation [iSj and func- 
tional renormalization based on an exact renormalization 
group equation 0, |23| for the effective average action 
pl| . We propose a partial bosonization which respects 
all symmetries. For a suitable choice the mean field the- 
ory (where only the fermionic fluctuations are taken into 
account) reproduces the Hartree-Fock result. For generic 
choices, however, the mean field phase diagram depends 
on unphysical parameters reflecting the "Fierz am- 
biguity" 22] in the formulation of mean field theory. This 
problem is overcome by the inclusion of the bosonic fluc- 
tuations. The size of the bosonic corrections depends in 
such a way on the choice of partial bosonization that 
the phase diagram becomes almost independent of the 
choice of the mean field already for the approximations 
employed in the present work |23l |. 



The scale dependence of the effective action is followed 
by the solution of a functional flow equation. This is for- 
mulated in the partially bosonized theory, employing the 
"rebosonization" technique of 24]. We thus combine the 
capability of mean field theory to describe an order pa- 
rameter with the advantages of the functional renormal- 
ization group. The renormalization flow includes step by 
step the bosonic fluctuations on larger and larger length 
scales and the fermionic fluctuations closer and closer to 
the Fermi surface as the "infrared cutoff" k is lowered. 
This formulation of the functional renormalization group 
can be used in the low temperature phase as welL in con- 
trast to the purely fermionic formulation 0, SS^S- 

A coherent picture for the antiferromagnetic phase at 
low temperature and small doping emerges. We find a 
simple effective action for the spin waves (cf. eq. (fT5|) 
and a more detailed discussion in which respects all 
symmetries. For temperatures near and below the effec- 
tive critical temperature Tc and for small doping (small 
chemical potential ^) it can be described by an 0(3)- 
symmetric effective potential U(a) for the antiferromag- 
netic order parameter a and an appropriate kinetic term 
which reads for small momenta {Zat^ /2)diddid. (Here Za 
is a wave function renormalization, t is the next neigh- 
bor coupling of the Hubbard model and the derivatives 
di are in lattice units.) For T < Tc the minimum of U 
occurs for ao = (Xq/2 7^ 0. The characteristic scale for the 
antiferromagnetic order is set by the renormalized order 
parameter Ka = Zat^ao/T . In particular, this determines 
the gap for the electrons, — 2T Ka / Z a{ha / 1) , with 
Yukawa coupling ha/t close to one for our numerical ex- 
ample 23]. Also the renormalized mass (or inverse cor- 
relation length) for the spin fluctuations in the "radial 
direction" (parallel to the order parameter ap) is deter- 
mined by Ka, ifia \/k^, whereas the spin fluctuations in 
the "Goldstone direction" perpendicular to ao are mass- 
less (infinite correlation length) for all T < Tc- Our result 
for the order parameter Kq is plotted in fig. 1. 

Actually, as a particularity of the two dimensional 
model, the size of the gap (and the value of Ka) de- 
pends on the infrared scale k which may be associated 
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FIG. 1: Renormalized expectation value Ka of the antifer- 
romagnetic condensate in the low temperature phase for 
M = 0, U/t = 3. 

to the inverse of the size of the experimental probe. (We 
take K, \cm.) This is due to the fluctuations of the 
massless Goldstone bosons which reduce Aq as k is low- 
ered. As we can see in fig. 2 this decrease is only logarith- 
mic. For any fixed T 7^ and ^ one always encounters 
a scale fc > where Kq vanishes. Therefore, the order 
parameter vanishes for fc — > 0, in agreement with the 
Mermin- Wagner theorem. Nevertheless, for low enough T 
the effective order parameter remains nonzero for any re- 
alistic macroscopic length k^^ . In this sense the Mermin- 
Wagner theorem fails to be practically applicable to the 
low temperature behavior of the two dimensional Hub- 
bard model. For fixed macroscopic k~^ and T < Tc the 
latter is well described by an ordered antiferromagnetic 
phase. 

The effective critical temperature Tc depends mildly 
on the size of the probe k~^ - a change of fc by a factor 
ten shifts Tc by 10%. Nevertheless, the two dimensional 
model exhibits effective critical behavior near Tc, with a 
macroscopic correlation length. For k <^ T the fermion 
fluctuations effectively decouple and the flow is described 
by a classical 0(3) linear scalar model (right part of 
fig. 121). For the correlation length f we find for T > Tc 

C(T) = ^exp|20.7/3(r)^| (1) 

with Tc = O.lUt , f3{Tc) = 1 and c{T),(3{T) smoothly 
varying functions of T ,23] without particular features 
for T 0. This closely resembles the result of '2^ which 
would correspond to constant /3. Near Tc the correlation 
length reaches the macroscopic size of the probe. Eq. 
is no longer applicable for T < Tc. 

The antiferromagnetic order occurs in the low temper- 
ature phase for a whole range of small doping. The critical 
temperature decreases mildly with the chemical potential 
fi 23] . We emphasize that the quantitative impact of the 
Goldstone-boson fluctuations on the value of Tc is very 
substantial. One may define a "pseudocritical tempera- 
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FIG. 2: Flow of the renormalized minimum of the potential 
Ka (solid), the quartic bosonic coupling 10~^Aa (long dashes) 
and the wave function renormalization 10^^ Za (short dashes). 
We have chosen U/t = 3 and T/t = 0.15; fi ^ 0. 

ture" Tpc where the mass term for the antiferromagnetic 
spin wave reaches zero for some fc, corresponding in the 
fermionic language to the four fermion coupling in the 
associated momentum channel growing "infinitely" large. 
We find Tpc/t = 0.18(0.21) for n = 0,Tpc/t = 0.145(0.18) 
for /i/t = 0.15 (with brackets for mean field theory) - to 
be confronted with Tc = 0.115 for fj, = 0. 

We next describe our formalism and approximations 
- for more details see We work within a functional 
integral formalism where the coordinates X denote the 
lattice sites as well as a "Euclidean time" t. We concen- 
trate in this note on fermion bilinears corresponding to 
spin density waves rn and charge density waves p. 

rh{X) ^ i,\X)ai,{X) , ~p{X) - ^t(x)^(X). (2) 

With these bilinears the four fermion coupling of the 
Hubbard model can be written as 

{i;\X)^{X)f = p{Xf = -\k{X)\ (3) 

We define the partially bosonized partition function by 

Z[ri,ri* ,Jp,Jm] = j V{ip*,'ilj,p,rh)exp{-{S + Sjj + Sj)) 

(4) 

with bosonic fields m, p, sources rj, J and classical action 

S S'_F,kin + ^Upp'^ + \U„i'm^ - Uppp - Ummm, 
Sr, ^ -r]^-tjj ~ ri'^ip*, S,j = -Jpp- Jm-m. 

In particular, the fermion kinetic term in momentum 
space takes the form {lof = (2n + l)7rT) 

5^^,kin = 5]V'^(0)Pi=^(Q)V'(Q), (6) 

Q 

Pf{Q) = iivp + ^ — f^, e{q) ~ ~2t{cos qx + cos qy) 

such that the Fermi surface is located at — 2i(cos(ji -t- 
008(72) = Here the Fourier transforms are defined as 
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V'(X) = ip{Q)e'^^^ and we use the short hand nota- 
tion Q = (t^^F, 9), X = {t,x), QX = ujpT + xq. 

The Gaussian functional integral over the spin density 
boson rh{x) and charge density boson p{x) can be per- 
formed explicitely. This demonstrates that Vl^ = InZ is 
equivalent to the free energy of the Hubbard model with 
next neighbor coupling and repulsive interaction U > 
(up to terms ~ which are explicitely calculable). The 
equivalence holds for all Up and U,n subject to the con- 
straint 

U = -Up + 3[/„. (7) 

The remaining freedom in the choice of Um is related to 
the Fierz ambiguity. The Hartree-Fock or lowest order 
Schwinger Dyson approximation corresponds to a com- 
putation of the fcrmionic fluctuation determinant (ne- 
glecting the fluctuations of p, m) for the choice Up = 
Um = U/2. (See JJj for a generahzed setting). 

The renormalization flow is implemented by the intro- 
duction of an infrared cutoff k through a piece in the 
action quadratic in the fields x — (V'jPj^) = ("01 ^)- 
presence of this cutoff the generating functional for the 
connected Green functions depends on the scale k 

Wk[J] = In J Pxexp ( - {S[x] + ASk[x]) + Jx)), (8) 

ASk = J2{4'\Q)Rt{QWQ) + lbi-Q)RliQ)biQ)}- 

Q 

For the fermionic cutoff we choose 

RtiQ) = iLUFi^ - 1) = 2m{nF + i)(Tfc - T), (9) 

which has the effect of replacing the temperature T by 
some function Tk in the fermionic propagator. We specify 
this function by = T'^ + k* , such that for large k the 
flow mimics a change in temperature whereas for k <^ 
T the modification of the fermion propagator becomes 
ineffective. For the bosonic cutoff we take 

R\{Q) = Z{e-Q^)Q{e-Q^), (10) 

where = in an appropriate interval, see eq. 1)16(1 . 
We define the effective average action as (x = dW/ dJ) 

Tk[x]^Jx-Wk[J]-ASk[x]. (11) 

The dependence on the scale k obeys an exact flow equa- 
tion [11 

dkVkbc] = iSTr{9fei?,[ri') + (12) 

which has a physically intuitive "one loop form" [2ll| . 
Here the "supertrace" runs over field type, momentum 
and internal indices and has an additional minus sign for 
fermionic entries. We see that the IR cutoff Rk adds to 



the exact inverse propagator matrix (second func- 
tional derivative) and therefore regulates possible zero 
modes, removing potential infrared problems for any 
fc > 0. The effective average action interpolates between 
the classical action S ("initial condition" for fc ^ 00) and 
the full "quantum effective action" for fc 0. 

As we are mainly interested in antiferromagnetic be- 
havior we define the boson (H = (0, tt, tt)) 

d{Q) = m(g + H) , a{X) = {~lf {X)<7^{X), (13) 

with P = 1 for odd lattice sites and P = else. The an- 
tiferromagnetic order corresponds to a{X) independent 
of X. 

For an approximative solution of the functional differ- 
ential equation (|12|l we employ a truncation of F^ con- 
sisting of a fermionic kinetic term, a Yukawa interaction 
between the fermions and bosons and a bosonic effec- 
tive action for the a-boson. We describe here the choice 
Up — Q , Um = C//3 whereas an extended truncation in- 
cluding the fluctuations of the p-bosons can be found in 
[2^ . For the fermionic kinetic term we adopt the classical 
part unchanged F^^^ = SF,kin{^)- The Yukawa coupling 
/ia,fe is taken to be scale dependent 

TY^k['4>,r,d] = -ha.k J2 a{KW{Q)(Ti^{Q') (14) 

KQQ' 

x5{K -Q + Q' + 11) 

and the purely bosonic part is described by a kinetic term 
and an effective potential 

Ta.k[d] = i^a(-g)P,(Q)a(Q)-|-^C/fe(a(X)). (15) 

Q X 

Due to S0{3) spin symmetry the fc-dependent effective 
potential Uk only depends on the rotation invariant com- 
bination a{X) = ^(i{X)a{X). The boson kinetic term 
{lub = 2'KmT) 

Pa{Q) = Za,kQ^ = ZaA^l+t^lQf), (16) 

involves a scale dependent wave function renormalization 
Za.k (with the dimension of mass~^) and the function [q]^ 
is defined as [q]^ — Qx^ly ^o^' 1i ^ I^""' ^^'^ continued 
periodically otherwise. 

The flow equation for the effective potential Uk{ct) ob- 
tains by evaluating eq. H12|l for a homogeneous antiferro- 
magnetic order parameter, a{Q) — (iS{Q), a = ^a^. The 
contribution of the fermionic fluctuations can be found 
by replacing in the mean field theory result T — > T/j in the 
fermionic propagator and applying the formal derivative 
dk = {dTk/dk)d/dTk- The bosonic contribution is the 
same as for the 0(3) linear cr-model pfl l2ll|: 

dkU{a) = -2T r ^dk\ncoshy{a) (17) 

+i ln[Fa(Q) + M2(a) -t- RtiQ)] 
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Here the squared fermion mass (gap) — 2h^a enters 
the function 



(18) 



The boson mass terms Mf obtain from the second deriva- 
tives of Uk and are different for the "radial mode" (in the 
direction of a) and the "Goldstone modes" (perpendicu- 
lar to a). In the SSB regime we recognize two massless 
Goldstone bosons for a at the minimum of the potential, 
a = Qfo 7^ 0. 

The flow of Za is determined by the anomalous dimen- 
sion ria = -kdk lli-Za, 



r,a = hltk{dkTk) 



[_d d_ 



(19) 



tanh^^+tanh '('>+''^^~^ 



Tfc7_, (27r)2 (e{q)+e{q + le,))/t 



1=0 



For the running of the renormalized Yukawa coupling 
ha — {T / ZaY^'^t~'^ha wc obtaiu in the symmetric regime 
(ao = 0) a direct and a "rebosonized" contribution 



kdkhl 



Phi 



girb) 
Phi 



-'^^'vf^'^ E {Pa\Q)PF\Q)PF\Q + n)} 

Q 

Q 

Pf\-Q)[Pf\'^~Q)-Pf\Q)]} (20) 



with ml = dUk/da\a=o , Pf^Q) = (ii^F - M + 
e(q))/t, LUF = (2n+l)^Tfe, p^iQ) = Ql/k^ + ml , Ql = 
Q^Q{Q^-k^)+k^Q{P-Q^). In the SSB-regime (ao 7^ 0) 
the change of the Yukawa coupling is negligible as we have 
checked numerically. Our results obtain from a numeri- 
cal solution of the flow equations (I17II . (|19|l . H2U|I with a 
quartic approximation to Uk(a) [23j. The initial condi- 
tions for very large k are given by ha = Um = C//3, Za = 
0, Uk{a) = Uma. 

We believe that our truncation catches the most impor- 
tant features in the region of small doping. Our method 
gives a unified description which covers at once several in- 
teresting limiting cases: perturbation theory and the per- 
turbative renormalization flow of the four-fcrmion cou- 
pling (small U), Hartree-Fock {Za ~ 0) and the classi- 
cal spin wave description (fc <^ T). The universal long- 
distance features (e.g. eq. (Q) are quite robust and should 
apply to a large class of models, including realistic mate- 
rials: for decreasing T the correlation length increases so 
rapidly that it reaches the macroscopic size of the probe 
k^^ at Tc- For T < Tc the formal symmetry restora- 
tion would require the randomization of antiferromag- 
netic domains with size much larger than k^^ . This effect 
is absent in a real experiment and macroscopic order can 
therefore be observed. 

The analysis can be extended to include other bosonic 
channels, including the superconducting channel relevant 
for large doping , . This would also lead to a more 
complete treatment of the momentum dependence of the 
effective four fermion interaction. 



[1] 



[2] 

[3] 
[4] 
[5] 
[6] 
[7] 

[8] 
[9] 



[10] 
[11] 

[12] 
[13] 

[14] 



J. Hubbard, Proc. Roy. Soc, A276, 238, (1963); J. [15 
Kanamori, Prog. Theor. Phys., 30, 275, (1963); M. C. 
Gutzwiller, Phys. Rev. Lett., 10, 159, (1963) [16 

B. I. Shraiman, E. D. Siggia, Phys. Rev. Lett. 61, 467 
(1988) [17 
R. Shankar, Phys. Rev. Lett. 63, 203 (1989) [18 
D. R. Penn, Phys. Rev. 142, 350 (1966) 
A. A. Gomes, P. Lederer, J. Phys. (Paris) 38, 231 (1977) [19 
H. J. Schulz, Phys. Rev. Lett. 65, 2462 (1990) [20 
D. Zanchi and H. J. Schulz, Z.Phys. B103, 339 (1997); 
Europhys. Lett. 44, 235 (1998) [21 

C. J. Halboth, W. Metzner, Phys. Rev. B61, 7364, (2000) 

C. Honerkamp, M. Salmhofer, cond-mat/0105218 [22 

C. Honerkamp, M. Salmhofer, T.M. Rice, 

cond-mat/0204063 [23 

L Grote, E. Kording, F. Wegner, cond-mat/0106604| [24 
P. van Dongen, Phys. Rev. Lett. 67, 757 (1991); 
Phys. Rev. B50, 14016 (1994) [25 

D. J. Scalapino, Phys.Rep. 250, 329 (1995) 
N.D. Mermin, H. Wagner, Phys. Rev. Lett. 17, 1133, [26 
(1966) 

F.J. Dyson, Phys. Rev. 75, 1736, (1949); J. Schwinger, 
Proc. Nat. Acad. Sc., 37, 452, (1951) 



J.W. Negele, H. Orland, Quantum Many-Particle Sys- 
tems, Addison Wesley, Redwood City, 1988 
T. Baier, E. Bick, C. Wetterich, Phys. Rev., B62, 23, 

15471, (2000) cond-mat/0005218 

T. Baier, E. Bick and C. Wetterich, |cond-mat/0107390] 

J. Hubbard, Phys. Rev. Lett. 3, 77, (1959) ; R.L. 

Stratonovich, Sov. Phys. Dokl. 2, 416, (1958) 

C. Wetterich, Phys. Lett. B301, 90, (1993) 

J. Berges, N. Tetradis, C. Wetterich, Phys.Rept. 363, 223 

(2002) hep-ph/0005122 

C. Wetterich, Nucl. Phys., B352, 529, (1991); Z. Phys., 
C57, 451, (1993); C60, 461, (1993) 

J. Jackel, C. Wetterich, hep-ph/0207094j, to appear in 
Phys. Rev. D 

T. Baier, E. Bick, C. Wetterich, cond-mat/0309715 
H. Gies, C. Wetterich, Phys.Rev. D65, 065001 (2002) 
hep-th/0107221 

G.v.Gersdorff, C. Wetterich, Phys. Rev. B 64, 
054513(2001) hep-th/0008114 . 

S. Chakravarty, B. Halperin, D. Nelson, Phys. Rev. B39, 
2344 (1989) 



